' V 


X 5 S/~ i 




X-621-|&-44 

; ; ' PREPRINT 


my., '■- 

j" ^ ^ ^ -x:. 


-^s, ,;V .~V. 


MUM 


w^. ; ; ;cv| 

DIRECT, IN : 

Sfc< -sj . . . j'y/ ■ 


V *\ 7 '• V&V \ 


V ' V ‘.fVcX 


OF WAVE PROPAGATION 
I N THE N EUTRAL THE RMO SPHERE 


• ,r -x ... - ... V. -x r.^v.', v ;• 

■ - '■' s v x ' Vy . .- / - / X t V 

- ■?- t -x 


w:.:;v 


,>. <- '■ r • 
■:( '•.;.. -/ 


7 .;f£ 


■ V ■ '. --, \~x - _ 

\Y ... ! -V ■ ■ J . - : 

,Y> VKV' y'- 

- » ^ r ^ ^ ) ‘ 


u u. ^ 


lxl J J 


,o 


:< .■ , x 

/ v. X 

- r 


' ■>: " l‘ Cl" 

■ ,> X. 

\ ; V'.x' i : ' 


P; NEWTON 


->,0 ■ ' ■ k V ' 

' \ y ' " 


rx 

3 s I s 

^1 K! 


zo 9 moi axiidvj 


w> </> 




O U. 


DAVID T. PELZ 'kJyf~ 
IANS VOLLAND , >f> V 

V ’ • -f v d V ' v , ; - - - V 

;.Cy ' v > v ,.; ! 




Y -h Cl 


V' -V- ? ; . • 

,i ■ v- ')>--* •)' ? ' ' '■ * ' 

y ■ ■// •■ ■ ■■ 4- -."'x v\ 


" > . ;■ i <- • A 

>x V'-- ^ 

\ - 1 ' ...,i > y ....... 


FEBRUARY 1968 


! •>~V N-. 


ut O a J 


\v * 'v-< - / 


O = 


o O 


-a u «g „ 

QC D. XJ 2 “ 

0 . _ 'z. o ■? 


- J • y t) 1: ■> -y V-- vc 

■ ' • ' • /'-.. 5 .'• .: '•• - 


9 w 

£L LL 

a o 


— O to ,- /t . , 

X 3 l- >. ~fT' : 

1-0 — 

co ~ 5 .^r 


i 2 


A, 


l-,„ '/• . 


t> ■ 

ISK> 


fc ' >v'^Vx. ; 


GSFC 


GODDARD SPACE FLIGHT CENTER 

V: GREENBELT, MARYLAND 




o vXn'.v-;^;- ' ? v ■ y . 

..V Yr 1 X ). 


i-M 


V^ : v^ 


> vV 


1 


X-621- 68-44 


DIRECT, IN SITU MEASUREMENTS OF WAVE PROPAGATION 
IN THE NEUTRAL THERMOSPHERE 


by 

George P. Newton, David T. Pelz, and Hans Volland* 


Aeronomy Branch 
Goddard Space Flight Center 
Greenbelt, Maryland 


‘National Academy of Sciences - National Research Council Senior Research Associate, on leave 


from Bonn University. 



2 


DIRECT, IN SITU MEASUREMENTS OF WAVE PROPAGATION 
IN THE NEUTRAL THERMOSPHERE 
by 

George P. Newton, David T. Pelz, and Hans Volland 

ABSTRACT 

The density gauges on the Explorer 32 satellite have measured local vari- 
ations in the atmospheric density which confirm that waves propagate in the 
neutral atmosphere. The waves were observed in the northern hemisphere over 
the altitude range from 286 km (satellite perigee) to at least 510 km. The waves 
are most prevalent at the higher latitudes near the auroral zone (orbit incli- 
nation is 65°) and were observed most frequently in the late evening and early 
morning hours, but were not limited to these latitudes and times. The virtual 
vertical half-wavelengths of the waves increase with altitude from 1 km at 286 
km altitude to 70 km at 510 km altitude. Four integrally related wavelengths 
have been observed, a "fundamental" and the second, third, and fourth harmonics. 
The wave half- amplitudes range from the limit of detectability to the largest ob- 
served of 50% of the smoothed density profile value. These waves are identified 
as free internal gravity waves propagating predominately north— south, or south- 
north, with horizontal wavelengths between 130 and 520 km and periods between 
10 and 100 min. This interpretation is consistent with observations of large 
scale traveling ionospheric disturbances by Thome and Georges. 
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DIRECT, IN SITU MEASUREMENTS OF WAVE PROPAGATION 
IN THE NEUTRAL THERMOSPHERE 
by 

George P. Newton, David T. Pelz, and Hans Volland 
INTRODUCTION 

It has been suggested by Hines (1960) that traveling ionospheric disturbances 
are caused by waves in the neutral air within the thermosphere. Previously 
these waves were observed only indirectly by the response of the ionospheric 
F-region to the neutral air disturbance (Heisler, 1967). 

This paper presents the first direct observations of neutral wave 
structure in the thermosphere and describes in greater detail than an earlier 
paper (Newton, et al., 1967) the characteristics of the waves. These results 
were obtained by the atmospheric density experiment on the Explorer 32 
satellite. Additional atmospheric density information is available from the 
experiment and will be presented in a later report. 

INSTRUMENTATION 

Satellite 

The Explorer 32 satellite, a second generation aeronomy satellite, was 
launched on 25 May, 1966 into a 64.6° orbit with a perigee of 286 km and an 
initial apogee of 2700 km. The satellite basic structure was a vacuum- 
sealed sphere , approximately 88 cm in diameter. The system had an active 
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life of approximately 10 months during which 8304 data producing turn-on’s, 
each of four minutes duration, were accomplished. A magnetic control 
system maintained the satellite spin-axis usually within five degrees of 
normal to the orbit plane and the spin rate at 1/2 rps. An artists’ con- 
ception of the satellite is shown in Figure 1 . 

The scientific instrumentation included in the satellite was: 

1. Two neutral particle mass spectrometers, one mounted on the satellite 
equator and one mounted on the spin axis ; 

2. One ion mass spectrometer mounted on the satellite equator; 

3. Two electrostatic Langmuir probes mounted on the satellite equator; 

4. Three magnetron type density gauges , two gauges mounted on the 
equator and one mounted 55° off the equator. 

The Satellite System was designed to permit turn on of the experiments and 
transmission or recording of data for approximately four minute intervals. Most 
of the experimental data were obtained in real time by the STADAN ground sta- 
tions, although considerable data were obtained at "remote" locations using the 
satellite tape recorder. 

Density Gauges 

The atmospheric density experiment used three cold cathode magnetron 
gauges of two types to measure the density. One gauge type (NRC-528) was 
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identical to the cold cathode gauges flown on the Explorer 17 satellite. The 
second gauge type (GCA-R5) differed in materials and geometry. Both gauges 
are shown in cross section in Figure 2. A description of the principles of mag- 
netron gauge operation is given by Redhead (1959). Descriptions of the flight 
gauges and their laboratory responses have been reported (Newton, et al., 

1963, Pelz and Newton, 1967). One gauge of each type was mounted on the satel- 
lite equator and the off -equatorial gauge was a GCA-R5. Each of the three gauge 
systems had an individual high voltage power supply and an electrometer. The 
equatorial GCA-R5 gauge had a linear, range switching electrometer which pro- 
vided high resolution-current measurements (> 98% resolution) and the other two 
gauges had logarithmic electrometers. All electrometers were calibrated once 
each satellite turn-on. The gauges were calibrated, sealed and launched under 
vacuum and opened on command when the satellite was in orbit. All instrumen- 
tation operated as intended in orbit. 

RESULTS 

The technique of using density gauges to perform atmospheric density 
measurements has been given (Newton, et al., 1965, Schultz, et al., 1948, Havens, 
et al., 1952, Horowitz and Kleitman, 1953). Measured atmospheric density is 
proportional to the pressure change measured by the gauges during the rotation 
of the satellite about its spin axis. Thus the time resolution of the density 
measurements by one gauge was two seconds (the spin period) corresponding to 
a spatial resolution along the satellite orbit of approximately 16 km. The wave 
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structure characteristics which were measured during a satellite turn-on are 
largely independent of the absolute calibration of the gauges and are influenced 
primarily by the relative precision of the gauge measurements which is believed 
to be better than 2%. Figure 3 shows a typical "smooth” atmospheric density 
profile. The line in the figure is determined by the least squares fit to the data 
of an exponential density dependence on altitude. It should be pointed out that the 
satellite was simultaneously changing local time, latitude and altitude during the 
pass. The changes in these parameters during a pass depend upon the location 
of the satellite in its orbit, however, typical maximum observed changes in local 
time and latitude during a turn-on are approximately 2 hours and 20 degrees 
respectively. 

Figure 4 shows the maximum and minimum "pressures" versus time meas- 
ured by the equatorial GCA gauge each satellite spin for tum-on 924. The local 
variations in the maximum pressure indicate structure in the atmospheric density 
distribution. During this pass the satellite was increasing altitude as evidenced 
by the decreasing gauge maximum and minimum pressure. 

Figure 5 shows the same data as in Figure 4, displayed as atmospheric 
density versus altitude. The local density structure is quite apparent and was 
observed simultaneously by the three independent gauge systems . The line in 
the figure again represents a least squares fit to the data of an exponential de - 
pendence of density on altitude. Complex wave structure is visible in Figure 5 
and suggests that there might have been a superposition of waves of different 
frequencies occurring. 
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Figure 6 shows the atmospheric density profile measured on turn-on 1067 
at altitudes near 400 km. 

Figure 7 shows the density profile measured at perigee during tum-on 637. 
The data start at 291 km altitude and proceed through perigee and back to 296 km 
altitude. The wave structure is apparent on both the increasing and decreasing 
altitude portions of the pass . 

Figure 8 shows the apparent vertical half -wavelength versus altitude de- 
termined from the passes listed in Tables 1 and 2. The lines are discussed in 
the interpretation section. The apparent half -wavelength has been plotted at the 
mean altitude of the wave half-cycle under consideration. The data have local 
times between 2100 and 0500 hours. The circles represent the shortest apparent 
half-wavelengths discemable during any turn-on and are usually easily meas- 
ured. The triangles represent the largest apparent half -wavelengths discemable 
during a tum-on when complex wave structure is present. The accuracy of the 
smaller half-wavelength determination is believed to be generally ±20%. The 
larger apparent half -wavelengths are more subjective in their determinations 
since some visual smoothing is performed. The two sizes of half-wavelengths 
appear to be integrally related, that is the smaller half -wavelengths usually 
correspond to harmonics of the larger half -wavelengths. The data have sur- 
prisingly small scatter, indicating that only a limited range of horizontal wave- 
lengths were observed (see discussion section). 
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The geographic locations of the passes where local structure was observed 
are shown in Figure 9, where the subsatellite points corresponding to the start 
and end data used from each satellite turn-on have been plotted and connected 
with a straight line. For the indicated days when wave structure was observed 
all data analyzed for that day are shown. It appears that the occurrence of 
the local structure is limited to latitudes greater than 20° and that the phenomena 
occurs more frequently near the auroral zone . 

The half-amplitudes of the density structure variations have been measured 
with respect to a line representing an exponential dependence of density on al- 
titude which has been fitted to the data by the method of least squares. The half- 
amplitudes vary from the limit of detectability to the largest observed of 50%. 

The most commonly observed half- amplitudes are 10-20%. At this time there is 
no consistent pattern apparent between the wave half -amplitudes and altitude, 
since the wave half-amplitude sometimes appears to be attenuated as the altitude 
increases, although constant amplitude waves are observed most often. There 
is some indication that the density is greater than usual during times when the 
local structure is apparent. The wave-like density structure is observed both on 
days indicated as geomagnetically disturbed and undisturbed. However, the wave- 
like density structure is observed predominantly during the late evening and 
early morning hours. 
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INTERPRETATION 

Traveling ionospheric disturbances have been observed which propagate with 
oblique phase planes and a bending of the phase planes toward the vertical with 
increasing height (Thome, 1964; George, 1967). Georges (1967), from HF-doppler 
studies, also observed vertically upward moving waves with periods of the order 
of 3 min which he interprets as acoustic waves. Hines (1967) has interpreted 
Thornes ? s measurements as resulting from free internal gravity waves which 
are refracted with height by the influence of heat conductivity within the 
thermosphere. Thus, two wave modes are suggested that could propagate in the 
neutral thermosphere. 

In this section we present first the conclusions drawn from the observed 
wave amplitude behavior. Second, we discuss the effect of the satellite velocity 
on the observed wavelengths. Third, we present observations of the relative 
locations of the observed wave phase planes in space and time, at the satellite 
perigee. Fourth, we discuss the implications of the preceeding three discus- 
sions on the type of wave mode observed and reach the conclusion that the 
observed wave structure is caused by internal gravity waves. 

1) Amplitude. 

Our experimental observations show that the amplitudes of the waves are 
nearly constant with height. From the discussion of the amplitude behavior of 
gravity and acoustic waves in the Appendix we conclude that the waves must 
propagate upward. Heat conduction wave amplitudes are attenuated more 
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strongly with altitude than gravity or acoustic waves and hence heat conduction 
waves can be excluded as the origin of our observed wave like structure. Thus, only 
gravity or acoustic waves remain as possible causes of the observed structure. 

2) The satellite velocity influence. 

For a theoretical consideration we assume that these waves are harmonic 
waves with angular frequency a> extending horizontally to distances of several 
thousand km with constant horizontal wave number k x . Assuming the validity of 
ray optics the phase of one of these waves is given by 


<&(x, z, t) 




Real (k z ) dg 


(1) 


where x is the horizontal direction of wave propagation, z is the altitude, t the 

* 

time and Real (k z ) is the real part of the (generally complex) vertical wave 
number k z of the dissipative wave within the thermosphere in which heat con- 
duction plays a dominant role. ( x x , z z , t : ) is an arbitrary initial point in our co- 
ordinate system. For convenience we select this point to be at the perigee of 
the satellite orbit. The projection of the satellite path into the (x, z)-plane 
leads to the following equations of motion: 


( x_x i) - Vxft-ti) 

(z-zj = V z (z)(t-t I ) 


( 2 ) 
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where V x andV z (z) are the horizontal and vertical satellite velocity components 
respectively. The y-component of the satellite motion does not influence the ob- 
served wave structure in our idealized model. The horizontal component of the 
satellite velocity is nearly constant and does not exceed the total velocity V: 

V = const = V .= 8.3 km/sec 

X 

The vertical component can be approximated near perigee by the analytic function 


/( z ~ z i) 


V z = ^ sign(t-tj) km/sec 


(3) 


with 


A = 19.8 sec/km 172 


We are interested in the vertical distance between two extrema of the wave 
structure along the satellite path. We find this apparent vertical half wavelength, 
A.', from the condition 

Z 



Eliminating x and t in Equation 1 by using Equation 2 we obtain 


z + V„/4 — z i 


z -X' /4-z, 1 

V z( Z -^/ 4 )J 


r z+x -V 4 

Real (k z ) d£ 

^ z -A. z /4 


77 (4) 
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and V z from Equation 3. From Equation 4 it is possible to calculate the apparent 
half- wavelength as a function of altitude. Equation 4 shows that the real vertical 
half -wavelength, which is defined as 

f Z+X z / 4 

Real(k z )d£ = 77 (5) 

•'z— X. 

is modified by an effect due to the horizontal motion of the satellite (the first 
term in Equation 4). This horizontal motion effect can lead to a negative co {{ , if 
co < k x V x and therefore to an apparent reversal of the direction of the observed 
phase propagation. 

The influence of the satellite motion on the apparent half -wave length varia- 
tion with altitude can be considered by using Equation (4) and neglecting the 
vertical phase integral. For a given value of gj ff the vertical half -wavelength 
may be calculated. The solid lines in Figure 8 represent the results of these 
calculations for values of <^> eff of 0.1 (line 1), 0.2 (line 2) 0.3 (line 3) and 0.4 
(line 4) rad/sec. It is seen that the altitudinal variation of the apparent half- 
wavelength can, for these values of cu , be explained solely by the variation 
of the vertical component of the satellite velocity with altitude. In Figure 8 it 
appears that the fourth harmonic is the greatest frequency observed. Apparently 
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the observed half -wavelengths for altitudes above 400 km. in reality belong to the 
fundamental half -wavelength. These results imply that the contribution of the 
phase integral in Equation 4 must be small. We thus have 

0.1 < oj , < 0.4 

err 

3) Relative Phase Plane Locations. 

Figure 10 shows the satellite versus time for the same data of Figure 7 
(turn-on 637). The location of the observed points of maximum and minimum 
phase of the wave are also indicated. The lines indicating planes of constant 
phase will be discussed later. The numbers in Figures 7 and 10 identify the same 
extrema in the density structure. However, in Figure 10 the uncertainty in the 
atltitude-time location of the extrema phase points are estimated by the dark 
'bands on the satellite orbit. As can be seen in Figure 10, the observations ob- 
tained on the decreasing altitude portion of the orbit are approximately 1000 km 
distance from the observations obtained on the increasing altitude part of the 
orbit. It should be noticed that amplitude and time spacing of the waves are 
approximately constant over the large horizontal extent of the observed structure. 
The time interval of 16 seconds between adjacent extrema in Figure 10 gives an 
apparent angular frequency of approximately 0.2 rad/sec. 

4) Observed Wave Mode. 

Two wave modes are considered as possible causes of the observed struc- 
ture, acoustic and gravity waves. We will now consider each type separately. 
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a) Acoustic waves. 

Contrary to the case of gravity waves, the thermosphere is nearly isotropic 
with respect to acoustic waves. Acoustic waves originating in the lower atmos- 
phere are expected to propagate mainly vertically because obliquely propagating 
waves are attenuated more heavily than vertically propagating waves by the in- 
verse distance dependence of the wave amplitude. This is in agreement with the 
measurements of Georges (1967) who observed the horizontal extent of coherent 
acoustic waves to be only a few hundred kilometers. There exists no observa- 
tional evidence of horizontally traveling ionospheric disturbances in the acoustic 

— 2 

frequency range (oo > 10 rad/sec.). 

The results of Figure 10 , which shows a wave-like structure with horizontal 
extent of more than 1500km, contradict an acoustic type wave interpretation. 
Moreover, for the equal horizontal spacing of approximately 16 seconds between 
adjacent extrema in Figure 10, and assuming k x = 0, one obtains a frequency of 
2 x 10 ~ x sec -1 . This frequency gives a vertical half -wavelength of 10 km for an 
acoustic speed of 0.7 km/sec. It is impossible to describe the extrema locations 
shown in Figure 10 by phase planes of constant frequency equally spaced in time 
with such small vertical wavelengths. Further, for small wavelengths the phase 
integral of Equation 4 would dominate the first term and a large difference in the 
observed virtual, vertical half-wavelengths between decreasing and increasing 
altitude satellite paths would be expected but is not observed. For these reasons 
we exclude acoustic waves as the origin of our observed wave structure. 
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b) Gravity waves. 

For gravity waves at these altitudes we have oo « <u eff and hence the apparent 
frequency is due to the horizontal motion of the satellite: 


iff 


A. 


k V 

X X 


277 V 


CO 


eff 


For co et - 0.1 

e r i 

A = 62.8V 

X X 

Therefore maximum A x £ (62.8) (8.3) % 520 km. 

For co e{{ = 0.4 

maximum A x 3:1 130 km . 

The maximum values of the horizontal wavelengths correspond to the wave 
propagating in or opposite to the horizontal direction of motion of the satellite 
and hence predominentlysouth-to-north or north-to- south. Georges (1967) and 
Thome (1964) report the predominent propagation direction of the waves they 
observe as north to south. This scale of horizontal wavelength also explains the 
nearly constant amplitude of the waves at higher altitudes since a small attenua- 
tion factor is expected for these wavelengths near 500 km (Volland, 1968). 

The large horizontal extent of the waves evidenced in Figure 10 is then quite 
consistent with the expected behavior of gravity waves. However, because the 
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magnitudes of the expected vertical wavelengths are large compared to the co- 
ordinate scales of Figure 10, the planes of constant phase in Figure 10 are nearly 
vertical and would be equally spaced in time. Thus the slope of the phase planes 
in Figure 10 cannot be determined from the points of extreme phase, and it is 
not possible to determine the real frequency of the waves from our observa- 
tions. From measurements of traveling ionospheric disturbances a frequency 
range between 10 -3 <w<10 -2 sec -1 is expected. As seen previously, the mean 
horizontal time spacing of the adjacent extrema in Figure 10 is approximately 
16 sec which corresponds to the second harmonic of oo efi = 0.2 sec - l . The 
vertical lines in Figure 10 have been drawn equally spaced by 16 sec intervals 
and are in agreement with the data. 

Equation 4 was used to calculate the apparent vertical half-wavelengths in- 
cluding the phase integral contribution for a horizontal wavelength of 520 km. 

—3 

and an angular frequency of 10 " rad/sec (period approximately 105 min.). 
Equation A1 from the Appendix was used to determine the integrand of Equation 
4, and an <w eff of 0.1 rad/sec was used. The results of this calculation are 
shown in Figure 8 as the dashed lines marked 1' and 1". The line 1' was calcu- 
lated for the satellite ascending in altitude and the line 1" was calculated for the 
satellite descending in altitude. It is seen that the vertical variation of the wave 
number and the vertical direction of motion of the satellite have little effect on 
the results from Equation 4. This is true not only for this particular case, but 
for gravity waves generally with frequencies 10 ~ < o> < 10~ . Thus, the 
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expected characteristics of gravity waves are entirely consistent with our 
observations. 

The real vertical wavelength is defined by Equation 5. Using Equation A1 
of the Appendix with Equation 5 the vertical half -wavelength, \J 2, was calcu- 
lated. Also calculated were the angles between the vertical and the directions 
of phase and energy propagation direction, d ph and d en respectively, from 
equations A5 and A6 respectively, and the attenuation factor 8 from equation A4. 
The results of the calculations are shown in Table 3 for altitudes of 300 km. and 
500 km. and two values of X. x equal to 260 km. and 520 km. with 
co = 10~ 3 sec -1 
H = 40 km 
<^ g = 7 x 10 -3 sec -1 
C s = 0.7 km/sec 
7 = 1.5 

f 7 x 10 ~ 2 sec -1 at 300 km altitude 

< _ 4 _ x 

(^7 x 10 sec at 500 km altitude. 

The values of S indicate that the attenuation of the wave amplitude is smaller 
for the larger horizontal wave length as suggested by Equation A4. For this 
wave length the amplitude falls to e -1 after the wave has traveled vertically 
upward approximately 300 km. It is also evident from the data in Table 3 that 
the wavefront and direction of energy propagation change with altitude and be- 
come nearly vertical above 500 km. altitude. 
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DISCUSSION AND CONCLUSION 

The consistency in the apparent wavelengths observed at various altitudes 
indicates that the same phenomenon is being repeatedly observed on different 
days. With this assumption it is possible to explain the observations on the basis 
of waves propagating in the neutral thermosphere with horizontal wavelengths 
such that the apparent altitudinal variation of half -wavelength is due primarily 
to the variation of the vertical satellite velocity. It then appears that the waves 
are gravity waves propagating predominently horizontally in a north-south or 
south-to-north direction: 

This conclusion is consistent with the results reported by Thome (1964) and 
Georges (1967) if it is assumed that the wave phenomena reported here is the 
underlying cause of the ionospheric effects with similar characteristics. This 
result supercedes the earlier preliminary conclusion (Newton, et al., 1967) in 
which purely horizontally propagating waves were excluded because of; 1) under- 
estimation of the doppler effects due to satellite motion; 2) the obscuration of 
extrema in the density data at perigee by local time changes and 3) an over- 
estimation of the accuracy of the extrema locations in time and space. 

The small effective frequency spread observed may result from possible 
filtering action of the atmosphere in addition to any frequency characteristics of 
the unknown wave source. It is interesting to note with respect to a wave source 
however, that in addition to an apparent relationship between the frequency of 
observations of the wave structure and proximity to the auroral zone, there 
might also be a relationship with the polar frontal zone cyclogenesis. 
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The energy propagated by the observed waves appears to be very small 
relative to the EUV energy flux into the thermosphere from the sun. The small 
energy flux could perhaps account for the waves being observed almost exclusively 
in the late evening and early morning hours . If the energy flux of the wave is con- 
sidered to be constant, then the density perturbation of the wave might be larger 
in a less dense atmosphere. Hence, the diurnal density variation could reduce 
the density perturbation to a value below the experiment detection threshold 
during the day. Further, it is possible that energy transmission characteristics 
of the atmosphere could perhaps change diumally thus diurnally changing the 
ability of the waves to propagate upward from lower altitudes. 
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APPENDIX 

The dispersion formula of neutral air waves within an isothermal thermo- 
sphere with finite heat conduction according to Pitteway and Hines (1963) is 
(using the notation of Volland (1968)) 


q - 


± 



+ 



"h\ 2 + 2 i "h 

CO } CO 




(Al) 


with 


q 




s 



CO 

a 




_ 7g 

"h " 2q 


c 


s 



(velocity of sound) 


K 

C h = c "jj a (heat conduction velocity) 


c 

7 = — (ratio between the specific heat at constant pressure 

V 

and constant volume) 

R ideal gas constant 
M molecular weight 
T temperature 
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g gravitational acceleration 

p density of the neutral air , 

H scale height 

k coefficient of heat conductivity. 

The plus sign outside the square root in Equation 6 is for downgoing waves, 
the minus sign is for upgoing waves. The plus sign inside the square root gives 
heat conduction waves, the minus sign gives acoustic-gravity waves. 

Equation A1 can be written 

±q = a. - j /3, (/3 = 0) (A2) 

where a is the real part of q and ft is the imaginary part. The real part is re- 
lated to the phase velocity of the wave, V, by 

C 

V = s (A3) 

/a 2 + S 2 

The imaginary part is related to the attenuation factor, e 8 ( z “ z o ) , of the wave 
during its vertical propagation by 

S = + k/? + ^1 (A4) 

For upward propagating waves the minus sign is used and for downward propa- 
gating waves the plus sign is used. Thus,- it is seen that S is always positive for 
downward propagating waves and that these waves are always strongly damped 
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in distances comparable to or less than two scale heights. For upward propa- 
gating waves, § is the difference between two positive values and the waves will 
be amplified or attenuated depending on the specific value of this difference. 

For gravity waves, using ray optics, the angle between the vertical direc- 
tion and the direction of phase propagation, d , is given by 

ph 

tan ■ (A5) 

p a 

The angle between the vertical direction and the direction of energy propagation, 
t> en , is given by 


tan r? 


en 


Ba 

"BS 


(A6) 


For altitudes greater than or near 500 km and co - 10~ 3 rad/sec, and S > 10 we 
have 



and therefore for upgoing gravity waves 
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Thus, 


tan i? 


ph 


CO 


2 coj3 

,ife 

; I CO 


(A8) 


indicating the well known downward phase propagation of upward propagating 
gravity waves. It is also seen from Equation A8 that the direction of phase propa- 
gation approaches the horizontal with increasing altitude because decreases 
with height. The attenuation factor becomes larger as the horizontal wavelength 
decreases indicating that the smaller horizontal wavelengths are attenuated 
more than the larger horizontal wavelengths at these altitudes. 
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Table I gives a summary of the altitude and time (from an arbitrary zero) 
of the extrema for turn-ons analyzed which show density structure. Also given 
are estimates of the horizontal and vertical satellite velocity components at the 
start and end altitudes corresponding to the density structure. 

Table II gives a s umm ary of the apparent vertical half-wavelengths and their 
mean altitude along with other pertinent information about the turn-on and density 
structure. 

Table III gives a summary of calculated propagation parameters for gravity 
waves with horozontal wavelengths of 260 km and 520 km at 300 km and 500 km 
altitudes. 
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Table I 





Maximum 

Minimum 

Start and End Satellite 
Horizontal Velocity 

Start and End Satellite 
Vertical Velocity 

Turn-On 

Number 

Date 

Orbit 

Altitude 

(km) 

Time 

(sec) 

(Arbitrary 

zero) 

Altitude 

(km) 

Time 

(sec) 

(Arbitrary 

zero) 

North-South 
+ to North 
(km/sec) 

East- West 
+ to East 
(km/sec) 

+ Upward 
(km/ sec) 

343 

6/15/66 

258 

288.9 

596.5 

287.4 

575.5 

6.83 

4.75 

0.022 






290.4 

611.5 







292.6 

629.7 

294.6 

644.0 







297.5 

662.5 

300.6 

678.0 







305.6 

701.5 



6.21 

5.49 

0.223 

371 

6/16/66 

272 



287.7 

535.7 

7.00 

4.49 

-0.089 




290.6 

503.5 



7.04 

4.42 

-0.118 

634 

6/25/66 

381 

288.4 

731.7 



6.79 

4.80 

0.070 

637 

6/25/66 

382 

289.3 

556.0 

287.7 

575.5 

7.13 

4.27 

-0.125 




286.8 

580.9 









286.1 

651.5 

287.4 

677.5 







288.7 

694.5 

289.9 

706.5 







291.8 

722.5 

293.1 

731.5 

6.69 

4.93 

0.109 

813 

7/1/66 

455 

292.6 

402.5 

295.6 

431.7 

6.75 

4.84 

0.114 




299.1 

454.5 

302.2 

462.0 







307.0 

483.0 

313.8 

506.5 







329.5 

564.8 



6.08 

5.60 

0.312 

924 

7/5/66 

504 

294.5 

857.5 

298.6 

881.5 

6.76 

4.83 

0.128 




302.4 

900.0 

306.5 

917.3 







309.8 

929.2 

313.9 

943.5 







320.0 

963.7 

326.7 

984.0 







336.6 

1013.5 

350.5 

1047.5 

5.68 

5.97 

0.407 

1066 

7/9 '66 

554 

286.5 

1066.0 

285.3 

1088.0 

7.05 

4.41 . 

-0.005 






287.5 

1054.4 







288.6 

1044.0 

290.2 

1030.2 







292.2 

1015.5 



7.24 

4.07 

' -0.172 

1067 

7-9/66 

554 

348.4 

423.5 



5.90 

5.76 

0.381 






385.0 

502.0 

5.25 

6.31 

0.487 

1074 

7/9/66 

557 

477.8 

536.5 

417.2 

441.5 

4.50 

6.80 

•0.572 






507.5 

578.1 

2.74 

7.58 

0.710 

1098 

7/10/66 

566 

327.9 

482.0 

320 

458.5 

6.26 

5.41 

0.303 






357.0 

556.3 

5.78 

5.87 

0.408 

1109 

7 '10/66 

570 

454.8 

574.8 



3.96 

7.09 

0.626 






526.0 

676.5 

2.34 

7.69 

0.735 

1479 

7/21 66 

702 

284.6 

312.0 

286.8 

339.7 

7.10 

4.34 

0.013 




292.4 

381.5 











294.7 

395.1 







297.3 

409.3 

300.4 

423.8 







306.8 

450.1 

313.2 

474.0 

6.64 

4.96 

0.239 
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Table H 


Turn-On 

Number 

Dale 

Average 
GMT 
(hr. min) 

Average 
Local Time 
(hours) 

Average 

Long. 

(degrees) 

Average 

Lilt. 

(degrees) 

Mean 

Altitude 

(km) 

Shortest 

2 

(km) 

Longest 

2 

(km) 

(gm cm') 

H. 

(km) 

Vertical Direction 
of Satellite Motion 

Geomagnetic 

Activity 

Wave 

Half- Amplitude 

343 

6 

15 66 

09 04 

5.82 

-48.7 

47.2 

288 

1.6 


1.4 14 

36. 2 

Ascending 

• 









290 

1.5 


1.3 14 

36.3 











292 

2.2 


1.25 14 

36.4 











294 

2. 1 


1. 2 14 

36.6 



10-15% 








296 

2.9 

8.9 

1.1 14 

36.7 











299 

3.6 


1.05 14 

36.9 











304 

4. 5 


9.5 14 

37.2 




371 

G 

16 G6 

12 04 

4.64 

-111.5 

31.8 

290 

3.4 


1.4 14 

36.3 

Descending 

_ 

10% 








292 


5.8 

1.3 14 

36.4 



634 

6 

25 66 

06 53 

3.60 

-49.3 

45,4 

288 


3.3 

2.2 14 

36.2 

Ascending 

D 

10% 

637 

6 

25 66 

08 47 

3. 30 

-83.0 

40.9 

288 

1.4 


2.1 14 

36.8 

Perigee Pass 

D 









287 

1.2 


2.2 14 

36. 1 











287 

1. 4 


2. 1 14 

36. 1 



10% 








288 

1. 4 


2.0 14 

36.2 










289 

1. 3 


1.9 14 

36.2 











291 

1.9 


1.7 14 

36.4 











292 

1.4 


l.o- 14 

36.4 




813 

7 

1 66 

05 59 

2.60 

-50.7 

48.0 

294 

3.0 


1.3 14 

36.6 

Ascending 










297 

3.5 


1.2 * 4 

36.8 











301 

305 

3. l 

4. 8 


1. 1 14 
9.7 •* 

37.0 

37.3 



10-20% 








310 

6. 8 


B. 14 

37.6 











316 


16.7 

7.4 14 

38.0 











321 

13.7 


6.4 

38.3 




924 

7 

5 66 

04 42 

1.76 

-44. 1 

48.0 

297 

4. 1 


1.15 14 

36.8 

Ascending 










300 

3.8 


1.05 14 

37.0 











304 

4.1 


9.0 14 

37.2 











308 

3.3 


8.7 

37.5 











312 

4. 1 


7.9 14 

37.7 



10-15% 








31G 

6. 1 


7. 1 " 

38.0 











323 

6.7 


G.0' 1 " 

38. 5 











332 

9.9 


4. 8 14 

39. 1 











344 

13.9 


3.5 ’’ 

39.9 




1066 

7 

9 66 

05 16 

23.7 

-83.5 

27. 2 

286 

1. 1 


2.3 14 

36. 1 

Descending 

D 









287 

1.0 


2.2 14 

36. 1 



10% 








288 

1.2 


2. 15 14 

36.2 










289 

1.7 


2. 1 14 

36.3 











291 

2.0 


2. 0" 14 

36.4 




1067 

7 

9 66 

05 24 

1.6 

-56.7 

54.4 

350 


27.5 

7. 0 14 

40.3 

Ascending 

D 

20-25% 







361 


24.0 

5.8 ,s 

41.0 










382 


37.0 

3.7 “ 

42.4 




1074 

7 

9 66 

11 15 

2.9 

-125.3 

61. 1 

472 

52 


S.O'"' 

48.2 

Ascending 

D 

20-30% 








511 

37 


2 1-16 

50.7 



1098 

7 

10 66 

04 36 

1.1 

-52.4 

51.7 

324 

8.2 


5.9' 4 

38.6 

Ascending 

D 

10-15% 








332 


18.0 

5. O' 14 

39. 1 










356 


30.0 

3. r' 4 

40.7 




1109 

7 

10 <66 

12 24 

3.2 

-137.6 

62.6 

528 

85 


1.7“ 

51.7 

Ascending 

D 

25-35% 

1479 

7 

21 66 

03 27 

21.9 

-83.2 

40.4 

286 

2. 1 


2.0' 14 

36. 1 

Ascending 

D 









294 

2.3 


1. 5' 14 

36.6 











296 

2.6 


1. 4‘ 14 

36.7 



10% 








299 

3.1 


1 . 3" 14 

36.9 










310 

6.5 


1.05 14 

37.6 
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FIGURE CAPTIONS 

Figure 1. Artist conception of Explorer 32 in orbit. 

Figure 2. Cross sections of the Explorer 32 magnetron gauges. The GCA-R5 
gauge is on top and has a stainless steel envelope. The NRC-528 gauge 
is on the bottom and has a glass envelope. Both orifice pieces are made of 
Kovar. No magnets are shown. 

Figure 3. Atmospheric density measured on turn-on 2762 over Rosman, North 
Carolina; on 25 August, 1966. The mean GMT of this pass was 19 hours 
34 minutes and the mean local solar time (LST) is 14.3 hours. The line 
represents a least squares fit to the data of an exponential dependence of 
density on altitude. The individual points are the density measurement made 
during each satellite spin by one gauge system. 

Figure 4. Maximum and minimum gauge pressures during each satellite spin 
measured on turn-on 924. The local variations in the maximum pressure 
indicate structure in the atmospheric density. 

Figure 5. Atmospheric density measured on turn-on 924 over Saint Johns, 
Newfoundland, on 5 July, 1966. The mean GMT is 4 hours 42 minutes and 
the mean LST is 1.8 hours. The line and points have the same description 
as in Figure 3. 

Figure 6. Atmospheric density measured on turn-on 1067 over Saint Johns, 
Newfoundland, at a mean GMT of 5 hours 24 minutes on 9 July, 1966. The 


mean LST was 1.6 hours 
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Figure 7. Atmospheric density measured at perigee during turn-on 637 over 
Rosman, North Carolina, at 8 hours 48 minutes GMT on 25 June, 1966. The 
mean LST was 3.2 hours. The numbers denote local maxima and minima 
in the density structure. 

Figure 8. Semi-logarithmic plot of apparent vertical half -wavelengths of waves 
propagating in the neutral atmosphere versus mean altitude of the half- 
wavelength. The solid lines marked 1, 2, 3, 4 are calculated from Equa- 
tion 4 neglecting the phase integral of Real (k z ) and using the values 
co e ff = 0.1, 0.2, 0.3, and 0.4 respectively. The dashed lines marked 1' and 1" were 
calculated from Equation 4 including the phase integral using co efi = 0.1 and 
a horizontal wavelength of 628 km and period of 6283 sec for the satellite 
ascending and descending respectively. 

Figure 9. Projections of satellite locations onto a world map for turn-ons for 
which atmospheric densities were measured on days when waves were ob- 
served. Turn-ons when waves were observed are dashed. The numbers 
beside the passes are orbit numbers. 

Figure 10. Locations on the satellite orbit of the maxima and minima of the 
density structure on turn-on 637. The numbers denote the same points as 
in Figure 7. The dark regions represent estimates of the possible spread 
in the observed location of the extrema of the density. The equally spaced 
lines of constant phase are drawn for a time separation of 16 sec. cor- 
responding to = 0.2. 
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Figure 1. Artist Conception of Explorer 32 in Orbit. 
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AUXILIARY 




Figure 2. Cross Sections of the Explorer 32 Magnetron Gauges. The GCA-R5 Gauge is on Top and Has a 
Stainless Steel Envelope. The NRC-528 Gauge is on the Bottom and Has a Glass Envelope. Both Orifice 
Pieces are Made of Kovar. No Magnets are Shown. 
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Figure 3. Atmospheric Density Measured on Turn-On 2762 Over Rosman, North Carolina; on 25 August, 
1966. The Mean GMT of this Pass was 19 Hours 34 Minutes and the Mean Local Solar Time (LST) is 
14.3 Hours. The Line Represents a Least Squares Fit to the Data of an Exponential Dependence of Den- 
sity on Altitude. The Individual Points are the Density Measurement Made During Each Satellite Spin 
by One Gauge System. 
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Figure 4. Maximum and Minimum Gauge Pressures During Each Satellite Spin Measured on Turn-On 
924. The Local Variations in the Maximum Pressure Indicate Structure in the Atmospheric Density. 
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Figure 5. Atmospheric Density Measured on Turn-On 924 Over Saint Johns, Newfoundland, on 5 July, 
1966. The Mean GMT is 4Hours 42 Minutes and the Mean LST is 1.8 Hours. The Line and Points 
Have the Same Description as in Figure 3. 
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ellite Ascending and Descending Respectively. 




— NO STRUCTURE * 6/16/66 o 6/25/66D * 7/5/66 ♦ 7/10/66D 

— STRUCTURE ° 6/15/66 ° 6/18/66 • 7/1/66 ■ 7/9/660* 7/21/66 D 


40 



36 


Figure 9. Projections of Satellite Locations Onto a World Map for Turn-Ons for Which Atmospheric Densities Were Measured on Days 
When Waves Were Observed. Turn-Ons When Waves Were Observed are Dashed. The Numbers Beside the Passes are Orbit Numbers. 











